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The results of a linear elasticity analysis yields that nano-rod inclusions aligned along the c 
axis of a thin filni of YBa2Cu307_5, such as BaZrOs and BaSnOs, squeeze that matrix by pure 
shear. The sensitivity of the superconducting critical temperature in that material to the latter 
implies that the phase boundary separating the nano-rod inclusion from the superconductor acts 



'^ ■ as a collective pinning center for the vortex lattice that appears in external magnetic field. A 

O . dominant contribution to the in- field critical current can result. The elasticity analysis also finds 



that the growth of nano-rod inclusions can be weakly metastable when the inclusion is softer than 
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Introduction. The ongoing development of thin films of superconducting YBa2Cu307_5 
(YBCO) for wire technology has resulted in world- record high critical currents.- In external 
magnetic field, the critical current is considerably enhanced by nano-rod inclusions that are 
aligned in parallel to the crystalline c axis.-'^'^^ The enhancement is strongest at orientations 
of the magnetic field parallel to the c axis. Understanding the fundamental physics behind 
this effect remains a challenge. It is also unknown what drives the growth of nanorods in 
the first place in YBCO films. 

In this paper, we provide insight into both of these questions by computing the strain 
field due to nano-rod inclusions that thread a YBCO superconductor along the c axis. The 
lattice constant of inclusions that optimize the critical current is typically 8% larger than 
that of the YBCO matrix in the a-b plane. Assuming a coherent phase boundary between 
the inclusion and a given epitaxial layer of YBCO, a linear elasticity analysis yields that the 
nanocolumn is compressed axially, while the YBCO matrix is squeezed by pure shear about 
the nanocolumn. The critical temperature in optimally doped YBCO is known to couple 
strongly to pure shear in the a-b plane.- In applied magnetic field, we show theoretically 
how this experimental fact results in substantial collective pinning of the vortex lattice by 
the phase boundary separating the nano-column inclusion from the YBCO matrix.-'-^ Also, 
the elastic energy shows weak metastability at a high density of nanocolumns when the 
nanocolumn is soft compared to the YBCO layer (see Fig. [1]). We believe that this drives 
epitaxial growth of nano-rod inclusions in YBCO films. 

Two- Dimensional Elasticity Theory. We shall determine first the elastic strain and the 
elastic energy cost due to a single nano-rod inclusion that threads a film of YBCO along 
the c axis. Such nanorods are typically composed either of BaZrOa (BZO)^*'^ or of BaSnOa 
(BSO).^ Both are cubic perovskites, with lattice constants (ain) that exceed that of the 
a-b plane in YBCO, flout = 3.86 A, by 9% and by 7% respectively.-^i^ Let us temporarily 
ignore the effect of the lattice mismatch along the c axis by considering only epitaxial 
layers that are far from any possible partial misfit dislocation, and that therefore present 
a coherent phase boundary between the inclusion and the YBCO matrix. Such partial 
misfit dislocations are accompanied by stacking faults,^ a topic which will be discussed later 
in the concluding section. The assumption of a coherent phase boundary is valid for a 
nano-rod inclusion of diameter less than the distance between possible misfit dislocations,- 
«Moire = («^ut ~ (^hi^)^^ ■ BZO uauorods typically have a diameter of- 2— 3nm, which satisfies 



the bound ctMoire = 5 nm. BSO nanocolumns, on the other hand, typically have a diameter 
ofi 7—8 nm. It exceeds OMoire = 6 nm, although not by much. 

Consider then a cylindrical nano-column inclusion that presents a coherent phase bound- 
ary with a given epitaxial layer of the YBCO matrix. Unit cells match up one-to-one across 
the phase boundary in such case. The ideal axial symmetry, assumed here for simplicity, 
implies a radial displacement field, u(r) = u{r)r. We then have the boundary condition 

Wout('"out) ~ '^in(^in) = ''"in ~ ''"out (Ij 

between the displacement fields of the nanocolumn (in) and of the YBCO layer (out) at the 
phase boundary. The in-plane lattice mismatch that it represents generates elastic strain in 
both the inclusion and in the YBCO matrix. The elastic energy due to a 2D strain field is 
given by the integral^ 



i^2D = /rfv|ic||(V-u)2 + icx (^ 



' dux duy\ ( dux duy^ 



(2) 



dx dy J \ dy dx J 

over the corresponding area (prime), which is confined to r < Tin for the nanocolumn and to 
^ > '"out for the YBCO matrix. Here, c\\ and c± are the 2D bulk compression modulus and 
the 2D shear modulus, respectively. A useful identity for the pure shear component above 
reads 

(|^-|fr-fe-y^-(-,^-,V.u,^-(Vxu,^. (3) 

The strain tensor takes the form Vu = {du/dr)rr+{u/r)(j)(j) in the present axially symmetric 
case. It combined with Eq. ([3]) results in the compact expression for the elastic energy, 
E2D = J d'^r{jC\\[r^^d{ru)/dr]'^ + ^c±[rd{r^^u)/dr]'^}. Calculus of variations then yields 
a nano-column inclusion squeezed by pure compression and a surrounding YBCO matrix 
squeezed by pure shear: 

M||(r) = —Aqt for r < rjn, and u±{r) = +BQrl^^/r for r > Tout, (4) 

with corresponding strain tensors 

Vu|| = -Aol and Vu^ = 5o(rout/r)^(00 - rf). (5) 

The total elastic energy ([2]) generated by the nano-column inclusion is then i?2D — 
2c\! nrf^Al + 2c °^ tct'^^^Bq. Minimizing it with respect to the constants Aq and Bq while 



enforcing the boundary condition ([T]) yields optimal values rin^o = {Ar)c"^ /(cy + c±^ ) 
and Tout -Bo = (Ar)cy /(cy + c\^ ). Here Ar = ri^ — rout- These then yield an elastic 
energy cost 

4^ = 2niArfic^f^-' + cf *^-')-i (6) 

for the nano-colunin inclusion, which has an equilibrium radius Tout + ""out ('"out) given by 

Consider next a field of many cylindrical nano-column inclusions of radius tq centered at 
transverse locations {Rn}- Suppose again that they all present a coherent phase boundary 
with a given epitaxial layer of the YBCO matrix. The displacement field is then a linear 
superposition of those generated by a single nano-column inclusion (jlj): 

Uin(r) = U|| (r - Ri) + ^ u^ [(rout/rm) (r - R,) + Ri - R^] (7) 

inside the i^^ nanocolumn, and 

Uout(r) = Eu^(r-R,) (8) 

j 

inside the YBCO matrix. The pure shear terms that have been added to the pure com- 
pression inside of a nanocolumn ([71) are required by the boundary condition ([T]). Ob- 
serve now, by Eq. ([5]), that V^uy = = V^u^. Inspection of the elastic energy 
functional ([2]) combined with the identity ([3]) then yields that the above superpositions 
are stationary because V ■ u_|_, V x uy and V x u_|_ all vanish. Indeed, the elastic en- 
ergy cost reduces to a sum of surface integrals around the phase boundaries of the form 
E2D = Ei -^20 + Ei Ej[ejj,i(out) + ei,ij(out)] + J^i Ej,fc[ei,j,fc(in) + eij,fc(out)], where the in- 
dices j and k refer to the terms in the superpositions (JTj) and (IHl), and where the index i 
refers to the phase boundary. The prime notation over the summation symbols indicates 
that i 7^ j, k. Each individual contribution Cij^k is given by a surface integral around the 
circle Si of radius rout that is centered at Rj: ejj_fe(X) = sgn(X) c]_ Ii,j,k) with 

/i,,- fc = i da- [Vui(r - R,)] • u^(r - R^). (9) 

^ Si 

Here, sgn(in) = +1 and sgn(out) = —1. Also, the measure da on the circle Si points radially 
outward. Substituting in the strain fields ([5]) above yields ultimately that ejj-.j = = ei^ij, 
and that 

e,,4X) = sgn(X)cf ) i27c)By^,,Re [i?,,,i?,,e^'^-'=« - r^]"' (10) 



for i 7^ j, k. (See Appendix.) Here, Rjj = Rj — Rj, and (pj,k{i) denotes the angle between the 
vectors Rjj and Rj,fc. The 2D elastic energy then is composed of a sum of 1-body , 2-body 
(j = k) and 3-body terms (j 7^ k), E2B = J2i -^'20 +Sj Y^'j,k ^i,j,ki with the interaction energy 
given by 

r,,,, = -{2'n)[c'r^ - cf\Blrl,Re [/2,,i?,,e^^^^'=» - rU~-\ (11) 

Notice that V^ j',a; changes sign as a function of the relative rigidity between the nano-column 
inclusion and the YBCO matrix. 

The elastic energy will now be obtained by computing subsequent self-energy corrections 
to the 2-body interaction and to the 1-body line tension. Let's first fix the coordinate 
for the phase boundary above, Rj, as well as one of the nanocolumn coordinates above, 
Rj. Observe that the 3-body interaction flTTl) has zero angle average about the center Rj 
over the remaining nanocolumn coordinate R^. This is due simply to the fact that the 
contour integral §dzz~^{z — w)^"^ around the unit circle, z = exp[i0j^fc(i)], vanishes for 
complex w inside of that circle. Let's assume that each nanocolumn has a hard core of 
radius ri ~ tq. At Rij ^ ri, we then obtain the estimate I]^ Vi,j,fc = — ^rr^ n<^Vi jj for 
the correction to the 2-body interaction on average over the bulk of the system. Here, 
Utf) denotes the density of nanocolumns. The renormalized 2-body interaction that results 
is then V^j = (1 — -Krln^)Vijj. Next, assume an effective hard-core of radius rg ~ ri 
for the nanocolumn at the coordinate Rj that remains. We thereby obtain the estimate 
Sj,fc ^j,fc = ^<f> S' d'^Ri,j^iJ = — 7rr2n0(l — 7rrJn(^)i?2D ^^ ^^e net self-energy correction to 
the elastic energy of an isolated nano-column inclusion, with 

rl = [(1 - cf VcrV(l + c^Vc^)] ■ [rtJir'i - r^)]. (12) 

This yields a total elastic energy density 

E2d/A = [1 - nrl n^{l - Txr\ n^)\n^E^S^ (13) 

as a function of the density of nanocolumns. The above third-order polynomial is depicted by 
Fig. [H It notably predicts weakly metastable epitaxial growth for relatively soft nanorods 
within the YBCO matrix, such that c\^ < c°^ . This occurs at a density n^ = (1 + 
[1 — {3rl/r2)Y^'^)/Snrf of nano-rod inclusions, at large effective crossections vrr^ > Sirrf. 
The equilibrium density of nano-rod inclusions therefore cannot be dilute. In particular, 
^nrln^ must lie somewhere between 1 and 2. Inspection of Eq. (IT2|) indicates that the 



former condition requires some degree of agglomeration among the nano-column inclusions: 
'"out < ^2 < 2rout- This may, however, be an artifact of the previous estimate for the 2-body 
self-energy correction, which is not accurate at Ri^ ~ ri. Last, the elastic energy cost per 
unit volume flT^ at meta-stable equilibrium is E2^/97irl = (2/9)(Ar/ri)^(C|| +0^^ )~^ 
in the marginally stable limit at r| = 3rf (see Fig. [T]). The strong dependence that it shows 
on the bulk compression modulus of the inclusion affects growth dynamics. This could be 
the root cause for the difference in length between BZO nanorods and BSO nanocolumns in 
YBCO.-^ 

Critical Current by Two-Dimensional Collective Pinning. We shall now determine the 
critical current of a thin film of superconducting YBCO threaded by nano-rod inclusions 
along the crystalline c axis and subject to external magnetic field aligned along the same 
axis. Recall that the critical temperature in an optimally doped YBCO superconductor 
is primarily sensitive to shear strain in the a-b plane.- That fact coupled with the shear 
strain generated by a nano-column inclusion ([5]) results in a potential-energy landscape 
for vortex lines that can collectively pin the vortex lattice. In particular, the contri- 
bution of the vortex core to the vortex line tension is approximated by the fundamen- 
tal energy scale per unit length Eq = ($o/4vrAL)^, where A^ denotes the London pene- 
tration depth. The temperature dependence shown by the vortex line tension is there- 
fore approximated by eo{T) = £:o(0)[l — (T/Tco)] near the mean-field critical temperature 
TcQ. The potential-energy landscape experienced by a vortex line then has a contribution 
fci(r) = J2aJ2i3{deo/dTc){dTc/dea,i3)ea,f3{T'), where Tc is the true critical temperature, and 
where ea,i3 is the symmetric strain tensor ([5]). It results in a ci-wave potential-energy land- 
scape about the nanocolumn for a vortex core, 

fei(r) = ep(rout/r)^cos2(/) , (14) 

with Ep = eo{0) (T/Tco)T~^[{dTc/debb) — {dTc/deaa)]Bo. Here the ratio between Tco and Tc 
is assumed to be constant. A rigid vortex line therefore experiences a force field 

fi(r) = /p(rout/r)=^(f cos20 + 0sin20) (15) 

due to the strain generated by a single nano-column inclusion, where /p = 2ep/ro^t is the 
maximum force per unit length. 

The above pinning/anti-pinning force (fT5|) is long range. The presence of an extended 
field of nanocolumns can cut the range off, however. (See Fig. [2l) Such forces add within the 



present elastic approximation ([8]): f(r) = J2i fi(r ~ R-j)- The d-wave nature of each isolated 
force field flTSj) implies a null net force on average. A characteristic fluctuation of the force 
over the YBCO matrix remains: /^ = n^ J' (Pr\fi{v)\'^ = ^(vrrQy^n0)/p, where integration 
(prime) is restricted to the YBCO matrix. Matching /^ with |fi(r)p yields an effective range 
for each pinning/anti-pinning center Vp = {2/7m^Y^^rJ^^^. 

The d-wave potential (TT^ that acts on rigid vortex lines in the vicinity of the phase 
boundaries between the nano-column inclusions and the YBCO matrix has zero angle aver- 
age. It therefore cannot pin down a vortex line in isolation. Previous work by one of the 
authors and Maley^ implies that many of them collectively pin the Abrikosov vortex lattice, 
however. A hexatic Bose glass state can exist at low temperature.- It is a vortex lattice 
threaded by isolated lines of edge dislocations in parallel to the relatively weak correlated 
pinning/anti-pinning centers. Plastic creep of the vortex lattice associated with glide by such 
edge dislocations limits the critical current,- which is given by jcB/c ~ fipf^/cQeb. Here rip 
denotes the density of vortex lines pinned by the nanocolumns, cqq = {^o/87i\L)'^nB is the 
elastic shear modulus of the pristine vortex lattice at a density ub of vortex lines,— and b de- 
notes the magnitude of the Burgers vector associated with the edge dislocations that thread 
the vortex lattice. The d-wave nature of the pinning/anti-pinning center flT^ also implies 
that its occupation is purely random. The density of vortex lines that they collectively pin 
is then equal to Up = {apnB)n^, where Up = 7r(rp — r^^J is the effective crossectional area 
of a pinning/anti-pinning center (see Fig. [2]). The critical current density therefore obeys 
a pure inverse-square-root power law with magnetic field, jc oc B^^^"^. Taking values of 
dTc/deaa = 230 K and dTc/debf, = —220 K for the strain derivatives of Tc in optimally-doped 
YBCO^ can result in a pinning efficiency, \fp\^/eo, of 93% at liquid nitrogen temperature! 

Discussion and Conclusions. We have found that the growth of nano-rod inclusions 
in YBCO films is very likely driven by weak metastability shown by the elastic energy of 
epitaxial layers. We also have pointed out how the sensitivity of the critical temperature in 
optimally-doped YBCO to pure shear strain inside of the a-b plane^ results in an effective 
collective pinning center for the Abrikosov vortex lattice at the phase boundary between the 
nano-rod inclusion and the YBCO matrix. 

The lattice mismatch along the c axis between the nano-rod inclusion and YBCO has 
so far been neglected, however. YBCO has a unit cell that can be divided into a stack of 
three cubes along the c axis, each with a lattice constant Cout/3 = 3.9 A. The strain that 



results at the phase boundary with a BZO nanorod or with a BSO nanocolumn, both of 
which are cubic with lattice constants ain = 4.2 A and 4.1 A, respectively, can be relieved by 
introducing partial misfit dislocations accompanied by stacking faults in the YBCO matrix.- 
The predicted spacing between such stacking faults, CMoirc = [(3/cout) ^(^h^]^^^ is then equal 
to 5nm for BZO nanorods and to 8nm for BSO nanocolumns (cf. ref. 5). Since their effect 
on the previous elasticity analysis can be accounted for by renormalized elastic moduli for 
the YBCO matrix, we believe that that our conclusions remain unchanged in their presence. 
The authors thank George Levin for discussions. This work was supported in part by the 
US Air Force Office of Scientific Research under grant no. FA9550-06- 1-0479. 

APPENDIX: SURFACE INTEGRALS 

Equation i^ gives the surface integral that determines the 3-body elastic interaction 
among nano-column inclusions. Integration by parts combined with V^u^ = yields that 
it is symmetric with respect to the latter: lij^k = h,k,j- In the case that ^ = j, it reduces to 
the angular integral 



1 /■2'r / R^, — r^ 



Ii.i.k = ^ri,,Bi j #1 '"^2/ -1|. (A-1) 



where rff^ = r^^^ + Rfi. + 2 rontRi,k cos(0 — (pk). Here, 0^ denotes the orientation of the vector 
Rj^fc = Rj — Rfc. After making the change of variables z = e**^, application of Cauchy's 
theorem yields that the integral vanishes: li^i^k = 0. In the case that i ^ j and i ^ k, the 
surface integral (jH]) reduces to 



T — _^4 r2 

^i,j,k — <'out-D0 ( 



2lT 



\^i,j ^out) 1 ^j,k _ \^i,j ^' out J Rj,k / ■ , x 

' i,j ' i,k ' ij i^k ' ij i,k 



■ (A.2) 



Repeating the previous steps results in a closed-form expression with a large number of 
terms. Symbolic manipulation programs then help reduce these to the result (fTOj) . 
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FIG. 1: Plotted is the total elastic energy density (fT3l) . in units of E'gD/^^i' versus the density of 
nano-column inclusions, in units of l/irrf. The dashed line above corresponds to the elastic energy 
of isolated nano-column inclusions. The radii in Eq. (J13p are set to r^ = =b3r^ for relatively soft 
and hard nano columns, respectively. 
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CHANGE IN VORTEX CORE ENERGY 
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FIG. 2: Displayed is a potential-energy landscape in units of |ep| and of the coherence length for a 
single vortex line that results from a superposition of 2744 d-wave collective-pinning centers [Eq. 
p^ ] arranged in a "liquid" fashion. (See ref. 8.) 
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